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Abstract 
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1^ We consider one-dimensional Pauli Hamiltonians in a bounded interval 

,__( with possibly non-self-adjoint Robin-type boundary conditions. We study 

,-H the influence of the spin-magnetic interaction on the interplay between 

(*7^ the type of boundary conditions and the spectrum. A special attention 

10 is paid to TT-symmetric boundary conditions with the physical choice of 

-Y~\ the time-reversal operator T. 
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1 Introduction 

In recent years there has been a growing interest in non-Hermitian "extensions" 
of quantum mechanics, usually associated with the names of yj- symmetry^ 
pseudo-Hermiticity, quasi- Hermiticity or crypto-Hermiticity (we respectively re- 
fer to [U [29l mj [38] where the first two works are recent surveys with many ref- 
erences). The quotation marks are used here because the extended theories are 
physically relevant only if the operators in question are similar to self-adjoint 
operators, which in turn puts the concept back to the conventional quantum 
mechanics. 

However, the freedom related to the existence of the similarity transforma- 
tion can be highly useful in applications, since a complicated non-local self- 
adjoint operator can be represented by a (possibly non-self-adjoint) differential 
operator (see |23| for one-dimensional examples), and the spectral theory for 
the latter is much more developed. Moreover, it is necessary that the non- 
Hermitian operators possess real spectra, which can be often ensured (at least 
in some perturbative regimes [^ I27|) by the simple criterion of CPT-symmetry. 

The goal of the present paper is to examine the role of spin in the above 
theories. We consider the simplest non-trivial situation of an electron (spin ^, 
mass m, charge — e < 0) interacting exclusively with an external homogeneous 
magnetic field B E R^. Choosing the Poincare gauge in which the magnetic 
vector potential coincides with ^B x x, this system is governed by the Pauli 
equation 

ih—- = A* + ^B-i*+ (Bx x)H + uB-a-9 =: H^ (1.1) 

at 2m h 8ni 

in the space-time variables {x,t), where h is the reduced Planck constant, fi :— 
he /{2m) is the Bohr magneton (for simplicity), L is the angular- momentum 
operator and c? is a three-component vector formed by the Pauli matrices. The 
spinorial wavefunction '^ can be represented as an element of L^(IR^) ® C^ and 



the operators appearing in ( 1.1 ) are assumed to appropriately act in this Hilbert 
space. 

The Hamiltonian H (equipped with a suitable domain) is Hcrmitian when 
considered in the full Hilbert space L'^{M.^^) (g) C^. Moreover, the Pauli equa- 



tion (1.1) is invariant under a simultaneous reversal of the space and time 
variables (c/ the discussion in Section pi). Relying on general definitions for 
the Dirac field (see, e.g., O §26]) and the fact that the Pauli equation can 
be obtained from the Dirac equation in a non-relativistic limit, the discrete 
symmetries can be represented by means of the parity CP and the time-reversal 
operator T (uniquely determined up to a phase factor). 



Our way how to "complexify" (1.1 1 is to restrict the space variables to a 



subset J7 C M and impose complex boundary conditions of the Robin type 

-—+A'$ = on dn, (1.2) 

on 

where n is the outward pointing normal unit to dil and A is a two-by-two 
complex-valued matrix. If Q is invariant with respect to the spatial inversion 7, 



it is possible to choose A in such a way that the yj-symmetry of ( 1.1 1 remains 
valid for the (possibly non-Hermitian) operator H on L^{^) (g) C~, subject to 



the boundary conditions (1.2) 



In this paper we study the interplay between the form of the matrix A and 
the spectrum of H. In particular, we are interested in the existence of real 
eigenvalues in the TT-symmetric situation. 

We are not aware of previous works on Pauli equation in the non-Hermitian 
extensions of quantum mechanics. However, there exist results on spinorial 
systems in the context of TT-symmetric coupled-channels models [351 IMl 133 
and the Dirac equation in the framework of Krein spaces [J [53] . 

One of the reasons for considering the spinorial model in this paper is the 
fact that the time-reversal operator T differs from the complex conjugation, the 
latter being the time-reversal operator for the scalar (i.e. spinless) Schrodinger 
equation, widely studied in the TT-symmetric quantum theory. In fact, for 
fermionic systems {i.e. half-integer non-zero spin), one has 

T2 = -1. (1.3) 

This has been remarked previously in the context of pseudo-Hermitian operators 
in [321 E] ■ A generalized concept of TT-symmetry as regards the operator 7 is 
suggested in [M) . 

The present model can be regarded as an extension of the one-dimensional 
scalar Hamiltonians with complex Robin boundary conditions studied in [211 
[2II1 E5] to the spinorial case. We refer to [32 [TS] for the discussion of relevance 
of (possibly non-Hermitian) Robin boundary conditions in physics and, in par- 
ticular, to Section [3] for a simple scattering-type interpretation in the present 
setting. 

This paper is organized as follows. In the following section we specify our 
model in terms of a one-dimensional Hamiltonian coming from ( |1.1[ ) . A physical 
relevance of the boundary conditions (1.2) is suggested in Section [3J Section [4] 



is devoted to a rigorous definition of our Hamiltonian as a closed operator as- 
sociated with a sectorial sesquilinear form. In Section [5] we discuss the physical 
choice of the operator TT and establish conditions on the boundary matrix A 
which guarantee various symmetry properties of the Hamiltonian. Section [6] is 
devoted to a spectral analysis supported by numerics; on several TT-symmetric 
examples we discuss the dependence of the spectrum on parameters character- 
izing the matrix A. The paper is concluded by Section [7] in which we mention 
some open problems. 

2 Our model 



We begin specifying our model represented by the Pauli equation (1.1 ). 

We choose the coordinate system in M.^ in such a way that the third coordi- 
nate axis is parallel with the homogeneous magnetic field B, i.e. B = {0,0, B) 
where B G M.. Then the orbital interaction B ■ L and the diamagnetic term 
{B X af)^ represent differential operators in the first two space variables only. 
On the other hand, the spinorial interaction B ■ a acts in the third space variable 
only (through the Pauli matrix a^ — diag(l, —1)). 

We set 

n := R^ X {-a, a) , (2.1) 



with some positive number a. Assuming that the matrix A in (1.2) is con 



stant on each of the connected components of dfl, the spectral problem for the 




Hamiltonian H therefore splits into two separate problems: a two-dimensional 
Landau-level problem in the first two variables and a one-dimensional problem 
in the third variable which we will study in sequel. Up to a constant factor repre- 
senting the energy of the given Landau level, the corresponding one-dimensional 
operators have the form 

Hb=\ dxy ^2 ,1 on J{:=i2((_a,a);C2), (2.2) 

~d?" 
subject to the boundary conditions 

*'(±a)+A±*(±a) =0. (2.3) 

Here we have put It? /{2m) — 1 and b := /^B, A^ E C^^^, and the third space 
variable is (with an abuse of notation) denoted by x. 

In view of the choice of physical constants made above, the only distin- 
guished length in our problem is the half-width a, and therefore the results 
must be scaled appropriately with respect to this length. In particular, the 
parameter b (characterizing the strength of the magnetic field) and eigenvalues 
of Hh (corresponding to quantum energies) become dimensionless when multi- 
plied by a?. The same can be done for the entries of A^ when multiplied by a. 
Consequently, all parameters can be thought as dimensionless in the sequel. 

As usual, the Hilbert space "K is identified with L^((— a, a)) (g) C^ and its 
elements are represented by the two-component spinors 

where V'± G ^^((^a, a)) (the ± notation should not be confused with the su- 
perscripts of the matrices A referring to the endpoints of (—a, a)). The inner 
product in 3-C is defined by 



r -T 
($,*) := / $ (x)*(a;)da;, 

J —a 



where the upper index T denotes transposition. The corresponding norm is 
denoted by || • || . The Euclidean norm of the spinor ^ as a vector in C^ is denoted 
by 1^1 :— -y/l'0+p -I- l'0-p and we use the same notation for the corresponding 
operator (matrix) norm \A\ := max {|AvI/| | ^1^ e C^, |*| = l} for A € C^^^. 

3 A scattering motivation 



Before giving a rigorous definition of our Hamiltonian formally introduced ( 2.2 ) 



(2.3), let us first justify the physical relevance of the boundary conditions (2.3 1. 

Our method is based on a generalization of an idea originally suggested in 15 1. 

Consider a generalized eigenvalue problem for the Hamiltonian of the form 



(2.2) on the whole space M locally perturbed by an electric field: 

^"'"'''^' . ° lfr)=Af*-). (3.1) 



Here x G M, A G M and V is the electric potential that is assumed to be 
compactly supported in {—a, a). Solutions ^ with A < — 1&| are bound states 
(associated with discrete eigenvalues), while those with A > —\b\ correspond to 
scattering states (associated with the essential spectrum). 

Outside the support ofV the problem (3.1 ) admits explicit solutions in terms 



of exponential functions. Let us look for special scattering solutions satisfying 



*(.t) = 



r^i^/X-fb X 



for 



Lt > a. 



(3.2) 



Then the (physical) problem (3.1) on the whole real axis can be solved by 
considering an (effective) boundary value problem in {~-a,a). The latter is 
simply obtained by considering (3.1) in {—a, a) and requiring that the solutions 
match at ±a smoothly with the asymptotic solutions (3.2). This leads to the 



boundary conditions (2.3) with an energy-dependent matrix 



At = 









-i^/ATfe 



(3.3) 



Note that (3.1) for x G {—a, a), subject to (2.3) with (3.3) at ±a, does not 
represent a standard spectral problem, it is rather an operator-pencil problem 
(because of the dependence of A^ on the spectral parameter A). It is non- 
linear in its nature. However, it can be solved by first considering a genuine 



(linear) spectral problem, namely (3.1) for x G {—a, a), subject to (2.3) with A^ 
at ±a, with a being treated as a real parameter. This leads to a discrete set 
of eigencurves a i->- A„(a), n G N. Then the "eigenvalues" of the true, energy- 
dependent problem are determined as those points A„(a*) satisfying the (non- 
linear) algebraic equations 

A„(a*) = a* . (3.4) 

The elements of the set {A„(Q;*)}„gN are called perfect-transmission energies 
(PTEs) in |15j , since their physical meaning is that they determine energies for 
which there is no reflection for the initial scattering problem (3.1) in E. It is 



interesting that PTEs are real, although they are obtained via solving a highly 
non-self-adjoint spectral problem. This feature is related to the fact that the 
choice A^ ensures that the boundary conditions are TT-symmetric (although 
not TT-symmetric in the context of the present paper where we do not allow 
the presence of A and b in the boundary conditions, see below). A physical 
interpretation of the possible complexification of the spectra of the auxiliar TT- 
symmetric spectral problem is also proposed in |15) . 

It is also interesting to note that switching on the static magnetic field 
{i.e. making b ^ 0) will typically lead to a splitting of the doubly degenerate 
eigencurves corresponding to the auxiliar non-self-adjoint spectral problem for 
6 = (c/ Figure l2|. Consequently, to each of the PTE in the scalar case 
without the magnetic field there correspond two PTEs in our spinorial model. 
The analogy with the Zeeman effect should not be surprising. 



The matrix (3.3) is complex and non-Hermitian, which is typical for effective 



models of scattering solutions of (3.1). On the other hand, real- valued Hermi- 
tian matrices are obtained when looking for bound states. In this paper, we 
proceed in a full generality by allowing arbitrary matrices A^ in (2.3). How- 



ever, it is important to stress that we regard the matrices as parameters entering 



the spectral system; the dependence of A^ on the spectral parameter A is not 
allowed and the dependence on the field b is allowed only if b is treated as a 
parameter (no change under the action of T, cf Section pi). 

To end up this motivation section, let us note that alternative proposals 
for the connection between non-Hermitian CPT-symmetric operators and physics 
have been suggested recently in the context of scattering in [THl [2H1 [3S1 E] ■ 

4 The Pauli Hamiltonian 

We now turn to a rigorous definition of the Hamiltonian formally introduced by 



(2.2)-(2.3|. In other words, since we are interested in spectral properties, we 
need a closed realization of the operator Hb- 

The easiest way is to define the Hamiltonian as the Friedrichs extension 



of the operator (2.2) initially considered on uniformly smooth spinors satisfy- 
ing (2.3). On such a restricted domain, an integration by parts easily leads to 
the associated sesquilinear form hb as a sum of three terms 

hb{<i>, *) ^ gi($, *) +bq2{^, *) + g3(*, *) , (4.1) 

where 

gi($,vl/) :=($', M/'), 

g2(*,*) :=($,CT3*), (4.2) 

g3($, *) := $^(a) A+ *(a) - #(-a) A' *(-a) . 

The form hb is well defined on a larger, Sobolev-type space 

D(M:=iji((-a,a);C2). (4.3) 

It is obvious for qi and q2, while the boundary term q^ can be shown bounded 
on 'D{hb) by means of the Sobolev embedding H^{{—a, a)) ^^ C'^d—a, a]). 

Our aim is to show that hb is a closed sectorial form. It is clear for qi defined 
on ( |4.3[ ), since qi is associated with the Neumann Laplacian {cf [101 Sec. 7]), 
and as such it is a densely defined, closed, symmetric, non-negative form. The 
term 52 represents just a bounded perturbation; indeed, |'Z2[^]| < II^IP for 
every ^P G JC. It is not longer true for q^, however, a suitable quantification of 
the Sobolev embedding can be used to ensure that q^ still represents a small 
perturbation in the following sense. 

Lemma 1. For every ^ e 'D{hb) and e G (0, 1), 

k3[*]| < e (i^+i + \A-\) \wf + C^^'+J^'' + '^^' + '^'' ) m' . 

Consequently, the form (72 + Qs *5 relatively bounded with respect to qi and the 
relative bound can be made arbitrarily small. 

Proof. The claim is based on the estimates 

|*(±a)|2 < 2||vl/'||||*|| + 1 ||vl,||2 < , |jvl,'|j2 +(^ + -) m' (4.4) 

Za \2a e J 

valid for any '^ e 23 (/ib). Here the first inequality can be established quite easily 
by the fundamental theorem of calculus and the Schwarz inequality. D 



Consequently, the perturbation result [T51 Tlim.VI.1.33] can be used to show 
that /if, is indeed sectorial and closed. According to the first representation 
theorem [181 Thni.VI.2.1], there exists a unique ni-sectorial operator H\j in 'K 
such that ft.6($,*) = ($,i?b*) for aU $ e ©(/i;,) and * S 2)(i/fc) C 1){hb). 
Following the arguments [IHl Ex. VI. 2. 16], it is easy to check that iJf, indeed 



acts as (2.2|-(2.3); more precisely 



Hu^ 






2) 



(i^b) = {* e H'^{{-a, a); C^) | *'(±a) + A±*(±a) == o} 



(4.5) 



Proposition 1. i/f, defined by (4.51 is on m-sectorial operator on "K. The 
adjoint of Hf, is given by 



Ht^ = 



-4,1 + byj+ 
-ij'f_ - biP^ 



'S{Hl) = [^^ e H^{{-a,a):,C') \ *'(±a) + (A±)* *(±a) = o} , 



(4.6) 



where A* = A^ . 

Proof. It remains to notice (c/ 18, Thm. VI. 2. 5]) that the adjoint operator 
is determined as the m-sectorial operator associated with the adjoint form hi 
defined by /i*($,*) := /i6(*, $), 2)(/i*) ;= S)(/ib). D 



Note that the choice A^ 
nian, subject to Neumann boundary conditions, that we denote by H^ . (At 
the same time, the choice A^ 
conditions.) 



gives rise to the (self-adjoint) Pauli Hamilto- 

idary conditions, that we denote by H^ . (At 

oo formally corresponds to Dirichlet boundary 



5 Symmetry properties 



It is well known that the Pauli equation ( 1.1 ) (in the whole space M ) is invariant 
under the simultaneous space inversion and time reversal (i.e. x i— >■ —x and 
t M- —t, respectively). This can be easily established if one realizes that the time 
reversal leads to a change of orientation of the magnetic field (i.e. B i-^ —B), 
while the orientation is unchanged by the space inversion. These properties can 



be deduced from Maxwell's equations to which the equation (1.1) is imphcitly 
coupled (c/ ESi §17])■ 
0ne is tempted to mathematically formalize the space-time reversal invari- 
ance in terms of a symmetry property of the Hamiltonian H. Given a unitary 
or antiunitary operator 6, we say that a linear operator _ff in a Hilbert space is 
C -symmetric if 

[i/,e] = o. (5.1) 

Here the commutator relation should be interpreted as an operator identity on 
the domain oi H, i.e. GH C HC In this framework, however, the Hamiltonian H 



appearing (1.1) is not TT-symmetric, just because there is no way how ensure 



the change of sign B under the action of T in the Hilbert-space setting (in 
which B is considered as an operator of multiplication). Nevertheless, H of 
course satisfies (5.1 ) with C = TT provided that the magnetic field is absent. 



One of the goals of this section is to determine the class of boundary matri- 
ces A^ which preserves the TT-symmetry in the sense above. In other words, 
since we do not like to think of 6 as a component of a field governed by the ad- 
ditional equations and to mathematically formalize the action of T on the field 
(6 is rather a fixed parameter in our Hilbert-space setting) , we restrict ourselves 
to rigorously looking for the property 

[Ho,y7]^0; (5.2) 



boundary conditions (2.3) satisfying this relation will be called TT-symmetric. 



In other words, boundary conditions arc TT-symmetric if, and only if, TT^ sat- 



isfies the same equations as '^ in (2.3). Similarly, we shall define TX-symmctric 
boundary conditions. 



In our one-dimensional situation (2.2), the parity T and the time reversal 



operator T act on spinors as follows (c/ [25, §30] and [551 §60], respectively) 



fp- (x) 



{y^)ix) := ^i-x) , (T*)(x) := ta2^ix) - ( ^^^ ) • (5-3) 

It is important to stress that T differs from the complex conjugation operator 



(3C^)(x) :=*(a;), (5.4) 

the latter being the time reversal operator in the scalar case. 

It is easily seen that T, T and X are norm-preserving, mutually commuting 
bijections on J{. CP is linear, while T and 3C are antilinear (i.e. conjugate-linear) 



operators. T and OC are involutive (i.e. T^ = 1 = 3C^), while T satisfies (1.3) 
Proposition 2. Hq is 

• "yj -symmetric if, and only if, A~ = 7A'^1, i.e., 

-^22 021 A r ^+ ^ fail ai2 

ai2 -all/ V"2i 022 

• "?%- symmetric if, and only if, A^ — —OCA^OC = — A+ , i.e., 

A-=(-^ -^) for A+=h' "1^ 

\-a21 -022/ va2i 022 

Proof. Since the space H^{{—a, a); C^) is left invariant under the actions of CP, T 
and 3C, it is enough to impose algebraic conditions on A so that the symmetry 
properties are ensured. More specifically, we need to ensure that ^ G !X)(i?o) 
implies TT^ £ 'I){Hq). Employing the identity 

(TT1')'(±a) + A±(J'T1')(±a) = (^7'^)'{Ta) + A^{7^)iTa) 

= -7 [^'(Ta) + 7A^7^{Ta)] 

and the bijectivity of T, the CPT-symmetry condition follows. The CP!}C-symmetry 
condition can be established in the same manner. D 



Another property we would like to examine in this section is related to the 
notion of S-self-adjointness. We say that a densely defined operator H on a 
Hilbcrt space is 5'-self-adjoint if 

H* = S^'^HS (5.5) 

for some bounded and boundedly invertible (possibly antilinear) operator S, 
where H* denotes the adjoint of H. It clearly generalizes the notion of self- 
adjointness and pseudo-Hermiticity. 

Proposition 3. Hq is 

• self-adjoint if, and only if, (A^)* = A^ ; 

• V- self- adjoint if and only if, A^ = — (A+)*, i.e., 

A- = ('^ '^] for A+ = f"" "12 

V-ai2 -a22/ \a2i 022 



yl='= = f " + ] with a^ e 



• 7 -self-adjoint if, and only if, [A^)* — — TA^T, i.e. 

3± 

a± 



• X-self-adjoint if, and only if, {A^)* = XA^X = A^ ; 

Proof. The claims follows by using similar arguments as in the proof of Propo- 
sition H D 

The spectral analysis of non-self-adjoint operators is more difficult than in 
the self-adjoint case, partly because the residual spectrum is in general not empty 
for the former. One of the goals of the present paper is to point out that the 
existence of this part of spectrum is always ruled out for 5-self-adjoint operators 
with antilinear S. 

Proposition 4 (General fact). Let H he a densely defined closed linear opera- 



tor on a Hilbert space satisfying (5.5) with a bounded and boundedly invertible 



antilinear operator S. Then the residual spectrum of H is empty. 

Proof. Since H is 5-self-adjoint, it is easy to see that A is an eigenvalue of H 
(with eigenfunction \1/) if, and only if, A is an eigenvalue of H* (with eigenfunc- 
tion S'""'^^). It is then clear from the general identity 

ar{H) = {A e C I A e cTpiH*) & A ^ crp(H)} 

that the residual spectrum of H must be empty. D 

The proposition generalizes the fact pointed out in [7] for S'-self-adjoint oper- 
ators with S being a conjugation operator {e.g. X) and applies to our (different) 
choice of T. 



6 Spectral analysis 

6.1 Location of the spectrum and pseudospectrum 

As a consequence of Proposition 111 we know that the numerical range of iff, is 
contained in a sector of the complex plane. Since the spectrum is a subset of the 
closure of the numerical range, it provides a basic information on the location 



of the spectrum of Hi,. However, coming back to the inequality (4.4 1 on which 
the proof of Lemma [T] is based, we are able to establish a better result in our 
case. 

Proposition 5. The spectrum of Hi, is enclosed in a parabola, 

a{Hb) C Sfe := jz e C I diz > - (\b\ + 4 \A\^ + \^\ =-. C , 

\Q-z\ < VS \A\ V^z + C + 1^ 

2a 

where \A\ := |yl+| + |A-|. 

Proof. By ^^ Corol. VI. 2. 3], the numerical range of Hi, is a dense subset of the 
numerical range of its form hi,, the latter being defined as the set of all complex 
numbers /ib[*] where * changes over all * G 2D(ft.b) such that \\'$\\ = 1. Using 
the first inequality of (|4.4[), we get 



^hi,[^] > gi[*] +6g2[*] - ksl*] 







> 


||^'||2 


- H 11*11' - 


2|A||1*'|1|1*||- 


1^1 
2a 


1*11' 






> 


^11-1 




2 

5 




9/l 


,[*]l 


< 


193 [*]| 


<2|A|||*' 


11*11+ ''^' 11*11^ 
2a " " 


) 





for every * £ 'I){hi,). The claim follows by combining these two estimates. D 

Thus the resolvent set of iff, contains the complement of Sf, in C As a further 
consequence, we can establish an upper bound on the norm of the resolvent: 

||(ii6-z)"^|| < l/dist(z,aSf,) for ah z e C \ Sf,. 

This result can be also interpreted as a location of the pseudospectrum of iff,, 
cf [III Sec. 9.3]. 

Remark 1. Note that the set Sf, in Proposition [5] is not symmetric with respect 
to the real axis. On the other hand, if iff, is C-symmetric with antiunitary C 
{e.g., if iff, is TOC-symmetric) , then we a priori know that the numerical range 
must be symmetric with respect to the real axis and an improved version of 
Proposition [s] holds. 
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6.2 The nature of the spectrum 

Since the Neumann Laplacian H^ has compact resolvent and the relative bound 
in Lemma [1] can be chosen less than 1/2 (in fact, arbitrarily small), it follows 
from [IHl Thm. VI. 3. 4] that Hi, has compact resolvent as well (for any choice 
of A±). 

Proposition 6. Hi, has a purely discrete spectrum (i.e. any point in the spec- 
trum is an isolated eigenvalue of finite algebraic multiplicity). 

Solving the eigenvalue problem Hi,^ = X^ consists in constructing the fun- 
damental system of —ip± — fc±'0± (say, in terms of sines and cosines), with 



k± :— \/X =F 6, and subject it to the boundary conditions (2.3). This leads to 
the following algebraic equation for the eigenvalues A: 

[det(y4+) + det(A^) — aiia22 — 022'^?!] ^-fc+ cos(afc„) cos(afc+) 
+ [det(yl+) det(y4^) + a^-^a'{-^k^_ + aj2«22^+ + ^-^+1 sin(afc_) sin(aA:+) 
+ [— det(A+)a22 + 0^2 det(y4^) + {—a'l-^ + a'^-^)k'^'\ k^ sin(afc_) cos(afc+) 
+ [— dct(A+)aj^j^ + a'^-^ dct{A^) + (—0^2 + '^22)^+] ^- cos(afc_) sin(afc+) 
+ («2iar2 + ai2a2"i) t^-K = 0, (6.1) 

where a^j and a~- denote the elements of the matrices A'^ and A~ , respectively. 



There are only a few choices of A'^ for which (6.1 ) admits explicit solutions. 
In the sequel we consider some particular situations that we analyse with help 
of numerical solutions. 

6.3 Examples 

A self-adjoint example with avoided crossings. Let us choose 

where a is a real parameter. It follows from from Proposition [3] that all the 
eigenvalues are real since Hi, is self-adjoint. The implicit equation for the eigen- 
values takes form 

2a^fc+fc_[l - cos(2afc+)cos(2afc_)] = -{k\kt + a'') sin(2afc+) sin(2afc_). 

The dependence of eigenvalues on the parameter a can be seen in Figure [l] An 
interesting phenomenon in this figure is the approaching of a pair of eigenvalues 
and its subsequent moving back and slowly approaching to constant values. It 
should be noted that in the point of closest approach the two curves do not 
intersect. This avoided crossing holds for each pair of the eigenvalues. 

A yj-symmetric example with real and complex spectra. As an ex- 
ample of non-Hermitian but TT-symmetric boundary conditions, let us consider 

,± (ia±P 



^ -1 0^ ialpl (6-3) 
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Figure 1: a-dependence of eigenvalues for b — 1 and a = ? in example (6.2 1, 



with a zoom of the avoided crossing of the first pair of eigenvalues on the right. 



where a and /3 are real parameters. The feature of this example is that the 
spinorial components do not mix. The implicit equation for the eigenvalues 
acquires the form 

{-2l3k^ cos(2a/c_) + {k'i ~ a^ - f) sin(2afc_)) 

X (-2/3fc+ cos(2aA:+) + {k\ - a^ - P^) sin(2afc+)) = 0. (6.4) 

Because of the decoupling, this eigenvalue problem can be analysed by using 
known results for this type of boundary conditions in the scalar case previ- 
ously studied in [21j and in more detail in [22] . It turns out that the spectrum 
significantly depends on the sign of f3. 

/3 = 0. It follows from [3T] that one pair of eigenvalues depend on the param- 
eter a quadratically and the others are constant, see the left part of Figure [2] 
More specifically, the eigenvalues explicitly read 



'^3,± - 



a^Tb if j=0, 

Tb if .?■ > 1 . 



jvr 
,1a 



(6.5) 



The crossings of full (respectively dashed) lines in the left part of Figure [2] cor- 
respond to eigenvalues of geometric multiplicity one and algebraic multiplicity 
two, while the crossings of full lines with dashed lines correspond to eigenvalues 
of both multiplicities equal to two. The entire spectrum is doubly degenerate 
for & = and there exist eigenvalues of geometric multiplicity two and algebraic 
multiplicity four. 

/? > 0. In this case, the reality of the spectrum was proved in [22 . The right part 
of Figure [2] shows the dependence of the eigenvalues on the parameter a. We 
again observe pairs of eigenvalues split because of the presence of the magnetic 
field. 

/? < 0. On the other hand, the reality of the spectrum in the case when /3 is 
negative is not guaranteed and, indeed, it is easily seen from Figure|3]that com- 
plex conjugate pairs of eigenvalues do appear when a couple of real eigenvalues 
collides as enlarging a. The pair of complex eigenvalues becomes real again for 
larger values of a. It follows from the analysis in |52] that only one pair of 
complex conjugate eigenvalues occurs simultaneously in the spectrum. 



12 




13 = 



13 = 0.5 



Figure 2: a-dependence of eigenvalues for b = 0.4 and a = j in example (6.3 1. 




Figure 3: a-dependence of the real (left) and imaginary (right) parts of eigen- 
values for b = 0.4, a = J and (3 = —0.5 in example (|6.3[). 



A TT-symmetric example with coupled spinorial components As an- 
other example of non-Hermitian CPT-symmetric boundary conditions, let us se- 
lect 

,± _ / ±ia^ 
~ l±ia 



A" 



(6.6) 



where a is a real parameter. The characteristic feature of this model is a non- 
trivial mixing of spinorial components. The implicit equation for the eigenvalues 
now takes the form 



4a^fc+fc_ cos(afc+)^ cos(aA;_)^ + 4a^fc+fc_ sin(a/c+)^ sin(aA:_)^ 

= — (fc_|_fc_ + a*) sin(2afc+) sin(2aA:_). 



(6.7) 



The dependence of low-lying eigenvalues on the parameter a can be seen in 
Figure |4] Here the lowest pair of real eigenvalues exhibits a crossing, however, 
the eigenvalues remain real after the crossing point as the parameter a increases. 
This behaviour is not featured uniquely by the lowest pair of eigenvalues, it also 
appears for higher- lying eigenvalues in the spectrum (not visible in the figure). 
On the other hand, as a increases, the other pairs of eigenvalues in the figure 
complexify after the first collision, then the corresponding eigenvalues propagate 
as complex conjugate pairs in the complex plane, meet again and become real. 
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ReX 





Figure 4: a-dependence of eigenvalues for b 
An animation can be found on the website [1 



0.5, a = v43 in example (6.6 1 



7 Conclusions 

The goal of this paper was to investigate the role of spin in complex exten- 
sions of quantum mechanics on a simple model of Pauli equation with complex 
Robin-type boundary conditions. A special attention was paid to CPT-synimetric 
situations with a physical choice of the time-reversal operator T. 

A simple physical interpretation of our model in terms of scattering was 
suggested in Section [3] It would be desirable to examine this motivation in more 
details and include "spin-dependent electric potential" {e.g. Bychkov-Rashba or 
Dreselhauss spin-orbit terms typical for semiconductor physics |14j). 

Robin boundary conditions represent a class of separated boundary con- 
ditions. Our model can be naturally extended to connected boundary condi- 
tions, whose spectral analysis represents a direction of potential future research 
(c/ mHHTg in the scalar case). 

In this paper we did not discuss the important question of the existence 

of similarity transformations (or the "metric" in the TT-symmetric context) 

connecting our non-Hermitian operators with self-adjoint Hamiltonians. The 

problem generally constitutes a difficult task and very few closed formulae are 

known (c/ [2(X, i3J [51 [23] and references therein). However, we can easily extend 

the results established in the scalar case without magnetic field [33] to our 

spinorial example (6.31 and compute the metric in this special case. Let us 

define 

^I + K 

I + K. 



9:= 



where / denotes the identity operator on i^((— a,a)) and K is an integral op- 
erator with kernel 

K{x, y) := ei"(^-!/)-/3|^-al [c + ia sgn(x - y)] , 
with c being any real number. It follows from [231 Sec. 4.5] and the nature 



of the decoupled boundary conditions (6.3 1 that & represents a one-parametric 
family of metrics for Hb under the yT-symmetric choice (6.3). More precisely, 
Hb is 0~^-self-adjoint (c/ (5.5)) and Q is positive provided that either: a is 
small; or /? is positive and large; or \c\ and |a| are small. To find the self-adjoint 
counterpart of Hi, determined by this similarity transformation constitutes an 
open problem (in the scalar case [23] there exists results for /3 = 0). 
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Our model was effectively one- dimensional. Higher dimensional general- 
izations in the spirit of [71 IS EO] would be especially interesting for variable 
boundary conditions (i.e. non-constant matrix A). 
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